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DEPART). ‘NT OF STATISTICS L) IC 


1. Introduction 


Cochran’s (1941) well-known test for equality of k normal population variances pro- 
ceeds as follows. Let o?,2 = 1,...,k, be the population variances, and the null hypothesis 
iS 

Ho : 07 = 02 =---=07 =o’, say. (1) 


Suppose k sample variances, one from each population. are given, each based on n degrees 
of freedom; let the values be s7,52,...,5%. Let these be ranked, so that, labelled in 
ascending order, the values are sit) < S{2) rn Six): To test Ho, Cochran suggested the 


statistic 


where Y = paw s?. Thus the statistic compares the largest sample variance with the 
sui of all variances, and is clearly appropriate to test Ho against the alternative that one 
population variance is larger than the rest. Solomon and Stephens (1989) have recently 
given very complete tables of the distribution of Z and have described, in particular, the 
usefulness of the test in process control; for example, to see if the largest variance in a 
sample of one week’s daily output is excessively large. Solomon and Stephens (1989) also 
give references to earlier literature and to other applications, particularly for n = 2, to 
testing for white noise in time series analysis. 

The question has been raised what to do if it is suspected that more than one vari- 
ance is particularly large compared to the others, perhaps, say, in examining 30 daily vari- 
ances from one month’s output. One might propose to apply Cochran's test successively: 
when one variance has been found too high, remove it and repeat the test on k — 1 vari- 
ances; continue this procedure, rejecting the largest of the current set of variances, till no 
further significance is found for Z. However, for such a technique, it is virtually impossi- 
ble to evaluate the overall significance level, that is, to find the probability of rejecting m 
variances, if this is what happens, when Hp is true. 

In this article extensive tables are given to test the hypothesis of k equal variances, 


each with n degrees of freedom, where rejection would indicate that the m-th largest sample 


variance is an outlier. This is based on the ratio of the m-th largest sample variance to ja35 


Special 


sVole ouu/Or 


the sum Y. Thus the test statistic is 


Steet ) 
Ms ee 
Z Y (3) 


Cochran’s statistic Z, defined in (2) is identical to Z,. The rationale behind this test is as 
follows. When the s? are ranked, we can imagine them plotted against the expected values 
of k ordered chi-square variables, with degrees of freedom n. If Ho is true, the plot should 
appear like a straight line. A commonly occurring alternative to the straight line plot, oc- 
curring in data analysis, is a plot in which the smallest s? values are close to a straight 
line, but, say, the top m values are all too large to be on the line. Suppose, then, the m-th 
largest s}?, that is s?,,,_,,): is significantly large. The other larger s?, namely s?,,5_ 1); 
Ste +3—m)> ete., will probably be large also, even if not themselves significantly large. be- 
cause they must be in order of increasing size. Such an event (s?, +1-m) Significantly large) 
will suggest that m population variances are larger than all the others. When, as usually 
happens, the value of o? in Ho is not known, it is necessary to divide by Y in order to 
obtain a scale-free statistic, and this leads to the statistic Z,,. The statistic is not formally 
developed as, say, a likelihood-ratio statistic, but will be useful in analyzing real data on 
variances, such as arise in quality assurance; in the example quoted above, several vari- 
ances during the course of a month may have appeared to be too large, and the test can 
be used to examine this possibility. 

The test procedure is set out below, and the thenry of the test, together with some 


comments, is given in Section 3. 


2. Test Procedure 


The steps in making the test are: 

(a) Rank the k given sample variances (which must be independent and each based 
on n degrees of freedom): s/,) < 8f.) < +++ < 8%). 

(b) Calculate Zm from (3), where Y = 7,) + 87.) +--+ 8?,). 

(c) Refer to Table 1 for the appropriate values of k, degrees of freedom n and m. 


(d) If Z is larger than the value given for significance level a. reject Hp at level a. 
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Table 1 contains the percentage points of the null distribution of Zm,. These have 


been found as described in the next section. 


EXAMPLE. Duncan (1986, p. 388) gives a data set of inside diameters, consisting of 
20 samples of size 5; hence k = 20 and n = 4. The values of s? are 16.5, 12.3, 10.3, 15.2, 
2.7, 11.3, 7.5, 19.8, 5.8, 5.8, 5.7, 11.3, 12.3, 3.5, 8.2, 12.5, 18.5, 6.5, 4.7, and 6.0; with a 
total of 196.4. This example was used in Solomon and Stephens (1989) to show the use of 
Cochran’s statistic (Z, in our present notation), and, with Z, = .1008, it was shown that 
the largest variance was not unusually large. It might be thought that the fourth largest 


sample variance 15.2 is an anomalous value, so we test the hypothesis of equal o? using 


Zs, which equals 15.2/196.4 = 0.0774. 


Reference to Table 1, with k = 20, n = 4 shows Z, to be not significant even at the 


50% level, so that there is no reason to reject Ho. 


3. Theory of the Test and Calculation of 


Percentage Points 


In Solomon and Stephens (1989) it was shown how the distribution of Z; could be 
calculated exactly for small k, and could be very well approximated, for larger values of 
k, by fitting Pearson curves. Again, for very small k, the exact distribution of Zm could, 
in principle, be found, although with greater difficulty, but here we use only the Pearson 


Curve approximation. Zm can be constructed as follows: 


(a) Let y1,y2,---, ye bei.i.d. random variables, each with the distribution o? ree where 


o2 


is any positive value; let y1) < y2) < --- < yxy be the order statistics of the set yi. 
(b) Let Y = Ujy;. 
(c) Then Zm = YW(k+1—m)/Y. 

It is clear that the distribution of Z is independent of a, the scale parameter of y,; also 


Y is a completely sufficient statistic for o?. Thus, by the Basu/Hogg/Craig Theorem, 


Zm and Y are independently distributed. We can henceforth assume that o = 1. Then 


4 


ZmY = Wk+1-m), and we have 


CNY pee 
2a) = Wier), (4) 


where E(-) denotes expectation. The denominator of (4) is easy to find, since Y is a y?- 


variable with kn degrees of freedom: then 


ep, 2° { (kn + 2r)/2} 


For the distribution of y(,41-m), let G(t) be the distribution, and g(t) the density, of y?; 
then E(yiy)s for any i (1 <i<k) is 


n! 


E(y(i)) = conics a t'G'""(t) [1 - G(t)]"~*g(t)de. (6) 


The moments of y(x41-~m) cam be calculated from (6), and hence, using (4) and (5), the 
moments of Zm can be found. The first 4 moments have been used to fit Pearson curves to 
the distribution of Z, as described by Solomon and Stephens (1978), and hence to obtain 
percentage points of Z,,. These can be expected to be very accurate, especially in the long 


upper tail which will be used, in general, in the present application. 


COMMENTS ON TABLE 1. (1) The case n = 2. For this case, the \3 distribution is 
essentially the exponential distribution, and the special properties of this distribution can 
be used to give an exact answer to the probability (Z, > z). This was done by Fisher 
(1929) in connection with testing for white noise in time series, and by Stevens (1939), 
in connection with a problem in geometric probability. Fisher (1940) discussed the two 
problems, and gave a small table of upper 5% points of Zm, for m = 1 and 2. (Z, and Zz 
are there called g,; and g2). A comparison of our points with Fisher’s is given in the small 
table following: P.C. refers to the Pearson curve points, F’ to the exact points as given by 


Fisher. 


k: 5 10 20 30 


m=1 F .6838*  .4449 .2704 .1978 
m= P.C. .6838*  .4437 .2700 .1977 
m=2 F .3670* = .2651 1755 .1336 
m= PC. .3653 .2641 1753 1335 


*These are exact results, see comment (2) below. 

It is clear that the Pearson curve fits give excellent results; calculations show that the 
errors in a given by using the P.C. points at the 5% level will be less than 0.25%. 

(2) For m = 1 an exact formuta is available for P(Zm > z) when z > 0.5. This was 
used in Solomon and Stephens (1989) and is again used in the present Table 1, which 
repeats results for m = 1 for the sake of completeness. 

(3) As a check, we calculated the moments of y441~m) using several available algo- 
rithms for the incomplete Gamma distribution. Table 1 as given uses the most recent, 
by Shea (1988). Other algorithms gave essentially the same values in the upper tail, but 
slightly different values in the lower tail, especially for small k and n. The lower tail is, of 
course, much less likely to be used in applications, and, to make the Table of manageable 
size, we give only the upper tail points. A much more extensive table including the lower 
tail, is available from the second author. 

This work was partially supported by the Natural Sciences and Engineering Research 
Council of Canada, and by the U. S. Office of Naval Research, and the authors express 
thanks to both of these agencies. 
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